Mesoscopic continuous and discrete channels for quantum information transfer 
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We study the possibility of realizing perfect quantum state transfer in mesoscopic devices. We 
discuss the case of the Fano- Anderson model extended to two impurities. For a channel with an 
infinite number of degrees of freedom, we obtain coherent behavior in the case of strong coupling 
or in weak coupling off-resonance. For a finite number of degrees of freedom, coherent behavior is 
associated to weak coupling and resonance conditions. 
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Quantum information and quantum communication 
require the ability of manipulating and transfer qubits 
in the space 0. Quantum state transfer (QST) can 
be realized by teleportation using flying qubits Q, 
or through quantum channels. When the information 
has to be processed in devices smaller than typical opti- 
cal "wavelengths for flying qubits, quantum channels are 
preferable. They can be based on solid state devices or 
on confined radiation fields. QST among optical cavities, 
as proposed by Cirac et al. some years ago |4{ , is possible 
due to the fact that each atom inside the cavity interacts 
only with a nearly monochromatic photon of the radia- 
tion field, and that photon can be transmitted unchanged 
to a distant site, before interacting with another atom in 
a second cavity. Generally, in mesoscopic devices, an in- 
teraction localized in the space involves all the modes of 
the support and the state reconstruction is affected by 
interference. 

The use of local excitations in quantum chains, first 
suggested by Bose Q, is indeed far from being optimal, 
due to quantum diffusion 0, 0] ■ Different physical real- 
izations of quantum channels have been suggested: fer- 
romagnetic spin chains 0, Q , Josephson arrays , na- 
noelectromechanical oscillators 9j. Quantum diffusion 
appears in each of these models. Ideally, this drawback 
can be overcome by using parallel chains and conditional 
gates 0] or through the adoption of en gine ered cou- 
plings between the nodes of the network |ll|. Practi- 
cally, these solutions are rather complicated to be real- 
ized. Recently, two schemes have been proposed in which 
the quantum chain act as a quantum bus, and both the 
encoding site and the decoding site are external and cou- 
pled locally with two different points of the chain [I3 , ll3j . 
In Ref . |l2| the coherent behavior of the channel is associ- 
ated with a simple model which takes into account only 
the center of mass mode of the chain. Numerical and 
analytical considerations show that the transmission in 
a finite-size chain can be nearly perfect and the speed 
of propagation is independent from the distance between 
sender and receiver, but only depends on the total length 
of the chain. 

The latter is one of possible models, where two (or 



more) quantum systems interact via an intermediate 
channel characterized by many degrees of freedom. 

In this Letter we study the general conditions under 
which such systems exhibit quantum oscillations, and 
therefore are suitable for QST, in the limit of discrete 
or continuous channel. 

For the sake of concreteness, we consider the Fano- 
Anderson model extended to two impurities: 



H = ^ £kc\c k + Q \ c\ca + c'gCB 
k 



(1) 



We have two quantum systems (A and B) with creation 
and annihilation operators c a ,ca,c b ,cb, a chain with N 
modes, described by ct (c&) which creates (annihilates) 
an excitation in the mode k, and interaction with the 
modes and A and B which amounts to tunneling pro- 
cesses in the case when both A and B are associated 
with a solid state tight binding model, or to a transfer 
of energy when A and B are atomic systems interacting 
with a radiation field. The coupling constant g mea- 
sures the strength of the interaction and the phase factor 
exp (ikL) takes into account the distance L between A 
and B. In the case of a continuous spectrum, sums must 
be thought as integrals. Due to the quadratic nature of 
the Hamiltonian, the evolution equation of each operator 
is independent from the corresponding quantum statis- 
tics. Then, the model works for fermions as well as for 
bosons. All the characteristics of the system are synthe- 
sized by the energy dispersion e^. 

In the case of continuum of states, possible candidates 
as mesoscopic channels are conductors in the tight bind- 
ing limit or one-dimensional wires with magnetic edge 
states 0], where there are experimental proofs of co- 
herent hopping with quantum dots 0, 0|. As far as 
discrete sets of states are considered, the model is suit- 
able to be implemented by arrays of quantum dots, or by 
nanoelectromechanical oscillators, or by a radiation con- 
fined in a finite-size cavity. An experimental evidence of 
coherent oscillations in an all solid state realization of a 
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J ay nes- Cummings- like scheme has been recently reported 
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We will show that coherent oscillations between A and 
B can be achieved using both continuous and discrete 
channels. In particular, discrete channels are suitable for 
our purposes when A and B are weakly coupled with 
the chain and Cl is resonant with one of its eigenvalues 
efc. In this situation, only the resonant modes play a 
significant role and the effective Hamiltonian is that of a 
few-body problem. The same behavior can be attained 



J 



with continuous channels in the case of strong coupling, 
or, in the weak coupling limit, whenever Cl lies outside 
the energy band. 

Let us start with an initial state where an excitation 
is present in the impurity A and both the second im- 
purity and the channel arc in their respective vacuum 
states: \ipi n ) = c\ |0). Writing the Heisenberg equations 
and replacing the operators with their Laplace transform 
c* (uj) = J °° e^'c^ (t) dt, assuming h = 1, we find 
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where 



A d (uj) 



N 
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and D(ui) = [ui-Cl- A (w)] 2 - A 2 L (ui). In Eq. ©, 
terms in c k , due to the presence of the excitation in the 
channel, introduce noise and limit the efficiency of the 
channel. 

Studying the zeroes of the spectral function D (uj), we 
extract all information about the system. We can assume 
ffe = — w cos ka, as usual when treating solids with many 
atoms and lattice constant a, and interpret A and B as 
impurity states. Here, 2w is the bandwidth and k is 
defined in the first Brillouin zone limited by and 2ir: 
k = 2itn/N , where n is any integer between and N — 1. 
Without loss of generality, we shall assume throughout 
the paper a — 1 and w = 1. It can be shown that, in this 
case, 



Ad (uj ) 



where Kd (uj) 
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K-d M + K N+d (uj) 



(uj 2 - 1) 



1/2 



1 



1 - K N (uj) 



1/2 



(4) 



In order to evaluate its zeroes, the spectral function can 
be decomposed in two factors: D(u>) = D + (uj)D_(uj), 
where 



D± (uj) = uj — f2 — 



1/2 



1 + K N (uj) ± [K_ L (uj) + K N+L (uj)] 
1 - K N (uj) 



(5) 



The analytic structure of Eq. (J2J consists in 2(N+1) real 
poles for every finite N, and, in the limit N — > oo, only 
4 real poles, related to the band extrema, remain, and 
poles inside the energy band are substituted by a cut. 



These exact results are compared, in the following, with 
the direct numerical solution of the evolution equation. 

We start from the weak coupling regime (g -C 1). The 
zeroes of Eq. (J5J can be calculated by iterating the zero 
order solution uj = ujq obtained in the limit g — > 0. 

First, we assume fi outside the energy band: |f2| > 1. 
In this case the zero order solution is ujq — Q, and, by 
iteration, 



wi, a = n + A (ft)±A L (O) 



(6) 



All roots are real and oscillations are expected. Residues 
associated to poles u>\ and ui% in Eq.(j2Jl are obtained ne- 
glecting terms in powers of order g 2 . In such limit we 
find that all the spectral weight is concentrated on the 
impurities' modes. Then, we obtain a coherent oscillation 
between the two impurities: 
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— c^-isin— c\ , (7) 



where uj + = 2 [CI + A (fl)] . 
In the limit of infinite 
lu+ = 2fl + 2g 2 / (VL 2 - 

2g 2 



and uj- 
number 

1/2 



= 2A z (n). 

of modes, 
1Y'~ and uj- = 

These solutions il- 
B experiences a Rabi 



[CI - Vfl 2 - 1J /y/SP - 1. 

lustrate that the open system A 
oscillation, and actually behaves as a closed one. Then, 
the system is suitable for QST or to create entangle- 
ment. In the case discussed above the dependence on 
the size-system is not crucial and the continuous limit 
is achieved even for not very large values of N. In Fig. 
H we report the probabilities of the excitation to be 
localized either on the first impurity or on the second 
one. 

The discussion becomes more interesting when |f2| < 1. 
In this case it can be useful to introduce an auxiliary 
complex variable 7 defined by u> — — cos 7, with the con- 
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FIG. 1: Time evolution of the occupation probabilities of A 
and B (Pa and Pb) in weak coupling and off-resonance. The 
coupling strength is g = 0.05, the impurities' energy is fl = 
1.5, the number of the channel's elements is N — 30, and the 
distance between A and B is L = 6. Here are reported both 
the numerical (exact) and theoretical curves. 



straint that < Re{j} < it. So 



D±( r y) = cos 7 — cosT 
1 
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sin 7 sin 7 ./V/ 2 



[cos 7 iV/2± cos 7(L + AT/2)], (8) 



having denned £1 = — cosT. At the resonance, fl coin- 
cides with one of the unperturbed poles. 

Since in the weak coupling limit the original energy 
levels are slightly modified, it is reasonable to assume 
that the resonant ones give the main contribution to the 
evolution and an expansion around them can be done. 
We write 7 = T + 5, with S expected to vanish in the 
limit of g = 0. Then 
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sin 2 T 



SN 

cot (1 ± cos FL) =F sin YL 



(9) 



Two different regimes appear for SN 1 or SN <C 1. 
In the first case the system is well approximated by 
its continuum limit, obtained replacing cotSN/2 with 
— i sign{Im{5}}. It is easy to show that Eq. 10 docs 
not provides polar solutions, but only singularities deriv- 
ing from the cut. Under these conditions, the excitation 
diffuses in the channel and the QST efficiency is lost. 

On the other hand, when SN -C 1 the 
cotangent in Eq.© is expanded into 2/ (SN) 
and sinTi is negligible. The solutions are 

then St 1 = 



±9 J 2(1 - cos FL) /N sin 2 F and 



Sf = ±<?y 2(1 + cos TL) /N sin T. The time evo- 
lution of ctt looks very simple when SI = and L is even: 



in such a case 
c\ (t) = cc 



-^4 + (-r/^/w^= 

fN A 



l+L/2 ■ 2 

' sin 



gt t 
l= c 

?N 



B 



+ k**%& +e -*)> (10) 



where ±fc are the modes in resonance with VL = 0. This 
formula shows that, despite the non vanishing probability 
of finding the excitation in the channel, perfect QST is 
achieved. In Fig. |2 we report the time evolution of Pa 
and Pb, which represent the occupation probabilities of 
A and B. On the other hand, assuming L odd, the second 
impurity is never populated: 



(t) = cosgt x / — c a+^ sin 9t\l — 



'-fe 



(11) 



The result of Eq. 1)11) |) is somewhat similar to that ob- 
tained in Ref. [12j, showing an efficiency of transfer in- 
dependent (limiting ourselves to even values of L) from 
the distance. The condition SN -C 1 (or gVN <C 1) can 
be interpreted as follows: the interaction splits the reso- 
nant pole in two levels with an energy separation of the 
order of g/y/~N, while the energy spacing between differ- 
ent modes is about 1/N. If none of the other modes falls 
inside this energy interval, then the excitation interacts 
effectively only with the resonant modes and the coher- 
ent behavior appears. Vice versa, when gj\[N ^> 1/N, 
the resonance is no more separated from the other modes 
and a continuum-like behavior is expected. 




FIG. 2: Numerical simulation of the evolution of Pa (t) and 
Pb (t) in weak coupling and resonance with the following pa- 
rameters: g = 0.01, Q = 0, N = 16, and L = 8. The time 
is normalized with respect to uj — \/2g/\fN. The theoretical 
behavior, calculated in the text, coincides perfectly with the 
numerical one. 

The other regime we want to explore is characterized 
by strong coupling (j> 1). Now, we look for solutions of 
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the order of g. then ft does not matter, at least in a first 
approximation, and can be set to zero. In practice, due 
to the strong character of the interaction, the modified 
energies of impurities lie ever outside the band, and the 
channel works in its continuum limit even for not too 
large N. Considering w > 1, we obtain, by iterative 
procedure: 



c A (t) — cos gt 



gt 



2(2ff) 



L A 



+ (-l) J 



i sin ■ 



gt 



2(2<?r 

dkf (k) 412) 



where / (k) is a function that satisfies the condition 
/ dk\f(k)\ 2 = sin 2 gt. In this case, we have high fre- 
quency oscillations between A and B and the channel 
modulated by a low frequency signal which enables QST. 
Note that the spectral weight is not entirely concen- 
trated on the impurities, because at intermediate times 
the probability of finding the excitation in the channel is 
finite. Note also the different scaling with the distance 
L of the information carrying oscillation with respect to 
weak coupling. In Fig. [3|the probabilities of finding the 
excitation on A and B are depicted. The lower panel 
shows the high frequency oscillation. The discussion of 
this limit fails when infinitely extended discrete spectra 
are considered, as, for instance, in the case of finite-length 
cavities. In this situation, coherent behavior is expected 
to come out only from resonance conditions. 

In conclusion, we have discussed by analytical and nu- 
merical calculations a number of possible configurations 
of a quantum bus allowing perfect state transfer or en- 
tanglement generation. The model we have considered 
is suitable for implementation of continuous or discrete 
quantum channels in different physical scenarios. We 
have illustrated in detail in which limits a coherent be- 
havior emerges, showing that both weak coupling and 
strong coupling are suitable for our purposes. We hope 
that our conclusions will help to design an experimental 
setup. 
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